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Abstract 

The oxidation of SiC-protected carbon/carbon through machined slots and naturally occurring craze cracks in the 
SiC was studied. The slot and crack geometries were characterized, and the subsurface oxidation of the 
carbon/carbon substrate at temperatures of 1000 to 1300 °C in air was assessed using weight change, x-ray 
computed tomography, and optical microscopy of sections. Rate constants were derived from these measurements 
and compared with a two-step diffusion control model of carbon oxidation. Oxidation kinetic measurements on both 
the specimens with machined slots and with naturally occurring craze cracks showed good agreement with the 
model. 
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1. Introduction 

1.1 Reinforced Carbon/Carbon 

Fig. 1 is a schematic of reinforced carbon/carbon (RCC) with a SiC conversion coating used on the Space 
Shuttle Orbiter nose cap and wing leading edges. The material is made with a two-dimensional layup of carbon- 
carbon fabric with repeated applications of a liquid carbon precursor to fill voids [1]. The resultant high strength and 
light weight make it an ideal aerospace material; however, oxidation is a major concern. An oxidation protection 
system is based on a SiC conversion coating. Because of the difference in coefficient of thermal expansion (CTE) of 
the SiC coating and carbon/carbon substrate, the SiC coating shrinks more than the underlying carbon/carbon on 
cooldown from the coating application temperature. This leads to vertical cracks in the coating, and these cracks are 
pathways for oxygen to reach the carbon/carbon substrate. Actual RCC used on the orbiters contains glass sealants 
to plug the cracks; however, this work examines carbon/carbon with only SiC. 

The crack morphologies in a similar system have been thoroughly characterized in an elegant study by 
Yurgartis et al. [2]. Both surface (plan) views and section views of the material were taken to obtain information 
such as crack spacing, crack width, and crack morphology. Such parameters were obtained for RCC in this study 
and prove critical to modeling oxidation. Buchanan and Little [3] have also characterized the crack pattern in a range 
of SiC coatings on carbon/carbon at various deposition conditions. 

1.2 Oxidation Through Coating Cracks 

Figs. 2(a-b) illustrate typical oxidation voids formed due to oxygen diffusion through a coating slot or crack, 
respectively, and subsequent reaction. A recent study by the authors outlines the microstructural features which 
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distinguish this void from a processing void [4]. In general, an oxidation void is associated with a path for oxygen 
ingress (i.e., the slot in Fig. 2(a) and crack in Fig. 2(b)). At higher temperatures (>1000 °C) these oxidation 
processes are expected to be diffusion controlled. Diffusion control implies a lack of preferential oxidative attack 
and hence uniform attack of the carbon/carbon with a plane interface at the void boundary, as illustrated in Fig. 1. 
However, a close examination of the voids in Figs. 2(a-b) indicates some preferential attack along fiber edges and 
the resultant “pointed fibers,” which are characteristic of reaction control [5]. Thus the actual process is a mix of 
diffusion and reaction control. It will be shown that diffusion control dominates. 

There are numerous studies of carbon oxidation through cracks and fissures in the protective coating and 
associated models found in the literature. Medford [6] developed models of oxidation for RCC, assuming pathways 
for oxygen are the cracks and fissures in SiC. His model encompasses the following processes: 

(1) Diffusion of oxygen inward through the cracks and/or fissures in the SiC 

(2) Oxidation of SiC crack walls to form Si0 2 

(3) Oxidation of carbon/carbon — matrix, fibers, or both 

(4) Diffusion of CO outward through the cracks/fissures in the SiC 

Medford developed total mass change expressions based on these processes. The expressions were the diffusion 
equations with modifications for the area of reaction. First he examined the case below 982 °C, where passive 
oxidation of the SiC is thought to be negligible; then he examined the case above 982 °C, where passive oxidation of 
the SiC is thought to be significant. 

The general process of carbon oxidation in a nonreactive matrix or through a nonreactive coating has been 
explored by several investigators [7]-[ 13]. Carbon oxidizes to CO(g): 

C(s) + 14 0 2 (g) = CO(g) (1) 

However, thermochemically the equilibrium between 0 2 , CO, and C0 2 favors the products in the following reaction: 

■/20 2 (g)+C0(g) = C0 2 (g) (2) 

Thus 0 2 (g) and CO(g) cannot coexist in similar amounts — they immediately will react to form C0 2 (g). However, 
the pairs of species in reaction (2), CO(g) with C0 2 (g) and 0 2 (g) with C0 2 (g), may coexist. Thus a two-step 
oxidation process is required to model carbon oxidation. The two reactions are 

1. At the carbon/gas interface: 


C(s) + C0 2 (g) = 2CO(g) (3) 

Note that C0 2 (g) is the oxidizing species. 

2. At a distance away from the carbon/gas interface: 

co(g) + y 2 o 2 (g) = co 2 (g) (4) 

In an earlier paper [13], the authors adapted this two-step diffusion-controlled oxidation model to cylindrical 
pinholes in the coating. Actual oxidation damage was assessed with weight loss measurements and area 
measurement of oxidation voids below the pinholes, using cross sections. At 1400 and 1000 °C hemispherical voids 
were observed on a macroscopic scale, indicating uniform attack and diffusion control. The two-step model 
described the process reasonably well at 1400 and 1000 °C. However, at 600 °C only minimal oxidation was 
measured, and reaction control dominated. Hence the model did not fit at this lower temperature. 

In the present study, the two-step oxidation model is further extended to describe a half-cylinder trough which 
forms under a slot or craze crack in the SiC coating. First a series of experiments is performed using machined slots 
in the coating. This allows testing the model with well-defined oxygen pathway geometry. Then the study is 
extended to actual through-thickness craze cracks in the coating. Oxidation damage is assessed with cross-sectional 
views, interrupted or real-time weight loss measurements, and x-ray computed tomography scans. The kinetics of 
carbon consumption is then compared with the model. 
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2. Experimental 

2.1 Characterization of Slots and Cracks in Test Specimens 

The tests described here were done on 1.91-cm-diameter RCC disks with a SiC coating on all sides. These disks 
were obtained from Lockheed-Martin Missiles and Fire Control in Dallas, Texas. Two sets of experiments were 
completed. Specimens with machined slots were used for the first set; specimens with naturally occurring craze 
cracks were used for the second set. Prior to oxidation, both types of specimens were characterized in order to 
quantify the area of exposed carbon. 

In the first set of specimens, artificial craze cracks with well-defined geometries were made with diamond 
blades (Keen Kut Products, Hayward, CA) of 0.25, 0.51, 0.76, and 1.02 mm thicknesses. These slots were cut to the 
SiC/carbon-carbon interface. These slotted specimens are shown in Fig. 3. One slot was cut in each specimen where 
the oxidation kinetics was followed with weight losses; however, two slots of different dimensions were cut in the 
specimens where the kinetics was followed using microstructure examination and/or x-ray computed tomography. 
Slot dimensions were carefully measured with a traveling microscope as well as directly from the cross sections and 
reported in the results section. Total coating thicknesses were measured from cross sections. Because of the 
variation in coating thickness, 20 random measurements were taken and averaged. 

The second set of specimens contained only the naturally occurring craze cracks. These were characterized with 
techniques similar to Yurgartis et al. [2]. Examination of the surface (plan view) did not reveal the complete crack 
pattern; however, grinding ~300 pm of SiC clearly revealed the crack pattern, which was then traced, as illustrated 
in Figs. 4(a) and (b). Image analysis software (Adobe Photoshop CS2 (Adobe Systems, Inc., San Jose, CA) with 
Fovea Pro 4.0 (Reindeer Graphics, Asheville, NC)) was used to determine the total crack length per unit area. The 
view of a section, as shown in Fig. 5, allowed measurements of crack width at room temperature. Several such 
sections were examined, and an average of 10 measurements is reported. Crack widths were measured near the 
surface and the small narrowing of the cracks closer to the carbon/carbon was ignored for this approximate 
characterization. The approximate area of exposed carbon/carbon per unit area of sample is determined by 
multiplying the crack width by the crack length per unit area. 

The measured values are listed in Table I and compared to the values reported by Yugartis et al. [2] and 
Buchanan and Little [3]. “Specimen A” of Yurgartis et al. is SiC-coated carbon/carbon with no matrix inhibitor; 
“Specimen B” is SiC-coated carbon/carbon with a boron-based matrix inhibitor. It should be noted that the materials 
examined by Yugartis et al. [2] and Buchanan and Little [3] are different from the RCC materials in this study. Here, 
the RCC had a thicker SiC coating applied by a pack cementation process than the other specimens in which the 
coating was applied by a chemical vapor deposition process. Yugartis et al. examined cracks that did and did not 
transverse the entire coating thickness. In this study only through-thickness cracks are examined. Differences in the 
specific CTE of carbon/carbon and the applied coating would be expected to have a major effect on the crack 
parameters. Thus the differences listed in Table 1 are not unexpected. It is interesting that despite the differences in 
crack width and spacing, the ratio of the exposed area of carbon to the surface area of the sample is roughly the same 
for all materials. 

2.2 Oxidation Treatments and Sample Analyses 

The samples were oxidized in air in either a box furnace or a vertical tube furnace. The kinetics of the box 
furnace exposures were followed with weight measurements every 0.5 hr for a total of 2.5 hr at 1200 °C. The craze 
crack specimens were suspended from a Cahn C-1000 (Cahn Instruments, Cerritos, CA) recording electrobalance 
by a sapphire fiber with a platinum wire basket into the vertical tube furnace. To avoid platinum silicide formation, 
the platinum wire was lined with alumina beads. Bottled air at 1 atm (101 325 Pa) was used at a flow rate of 
1.67x10^’ liter/s, and isothermal oxidation tests were performed at 1000, 1 100, 1200, and 1300 °C. 

Selected samples were mounted and sectioned for optical microscopy to obtain measurements of oxidation void 
size. A specialized sample preparation procedure involved a premount of the entire sample in epoxy. A standard 
vacuum infiltration was done followed by a high-pressure application to fill the open porosity. Then the sample was 
sectioned, remounted in a standard metallographic mold, and polished with diamond pastes. This ensured clear 
microstructures and preserved the interface between the SiC and the carbon/carbon. 

X-ray computed tomography (CT) was used to provide additional images of the oxidation damage, without 
sectioning. This SmartScan Model 100 (CITA Systems, Inc., Pueblo, CA) system utilizes a Feinfocus FXE-160 
(COMET AG, Flamatt, Switzerland) microfocus x-ray source to produce very high-resolution imaging of samples, 
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approaching 25 pm, in the CT mode of operation. The major hardware components of this system included the x-ray 
source, an area detector system, a five-axis object positioning subassembly, and a lead-lined radiation cabinet. A 
dual-processor computer system controlled the data acquisition and image processing. The slice plane thickness was 
120 pm per slice for these samples. Putting together slices electronically gave a three-dimensional view of oxidation 
damage. 

3. Diffusion-Control Oxidation Model 


3.1 Diffusion through Slot or Crack Alone 


Fig. 6 is a schematic of the two-step diffusion-controlled model for oxidation through a slot or crack. Note that 
CO 2 is the oxidizer and reaction (3) occurs at position x = 0 whereas reaction (4) occurs at position x = Xf. A list of 
symbols used in this report is in the appendix to aid the reader. The following derivation is an adaptation of this 
standard two-step diffusion controlled model [7-13] to the coating slot or crack geometry. 

The molar flux, J f ( in weight/unit area-unit time), of each species, i, is defined by 


J,=D 


eff 


r dc^ 

dx 


. ave 

+ V; C: 


(5) 


The first term is the diffusive term and the second is the convective term. Here D, eff is the effective gas phase 
diffusivity, c, is the molar concentration of species i (i = CO, C0 2 , or 0 2 ), and x is the distance down the slot or 
crack. The average molar velocity vf ve is defined as 


2 >, Yj j ‘ 

i i 

L 'T 


Here c T is the sum of concentrations. Equation (5) becomes 


Ji=Df 


f dc A 


\ dx j 


C T ; 


The boundary conditions are shown in Fig. 6(b): 

At v = 0 : c co = c ° 0 c c q 2 = c » 02 

At X = Xf\ Cco 2 =C* co Co 2 =c co =0 


At x = L: c 0 , = c ( L h c C o 2 = 0 


( 6 ) 


(7) 


( 8 ) 


As shown in Fig. 6(b) L is the total crack depth, which is the coating thickness. 

Consider the fluxes in Region 1 (Fig. 6(a)). Here and hence the convective term in equation (7) is 

zero. The fluxes of 0 2 (g) and C0 2 (g) in this region are thus calculated as 


J b 2 =- D < 


o 2 


r dcp 2 

dx 


~ D o A, 

(■ L ~ x f ) 


(9a) 
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(9b) 


J co, d co 2 


r dcco^ 
dx 


°C0 2 c co 2 

( L ~ x f) 


A determination of the CO 2 flux emerging from the crack or slot, , is the ultimate goal, because it is a measure 

of carbon consumption. First the quantity xyis needed. This can only be obtained as the ratio XfIL, which is derived 
from the following equations (eqs. (10) to (17)). The two fluxes in equation (9) are equal and opposite. From this 
fact, one can deduce the concentration gradients in Fig. 6(b). The concentration of C0 2 at x = Xj us 


'CO, 


D <h c o 2 

D , 


CO, 


(10) 


Equation (7) for CO and CO, in Region II (Fig. 6(a)) is more complex, as now there is a convection term since the 
stoichiometry of reaction (3) requires 


J 


11 . 

co - 


-2J 11 

J co 2 


(11) 


Thus equation (7) for CO, becomes 


^co, D[ 


co. 


5c ( 


CO, 


dx 


c CO, J ci 


CO, 


C T 


( 12 ) 


Rearranging to separate variables gives 


J C o 2 ^ x d co 2 


f \ 

C T 

C T + C C q 2 


dc, 


C0 2 


(13) 


jp n 

This can be integrated from cqq 2 = Cqq 2 atx=xj and cqq 1 = Cqq^ at x = 0 : 


/"o =^£2^1nl 


Ct+C . 


CO, 


Ct+C , 


0 

CO, 


(14) 


* 0 

This can be simplified further using equation (10) and the fact that c'co, » : 


J 11 
CO, 2 


d co 2 c t 

x f 


In 


D 0 2 c 0 2 
1 + — 

d co 2 c t 


(15) 


A mass balance between CO and O, gives 


Jco=-2J ] 02 


(16) 


Combining this with equation (11), 


J 11 =J l 

J CO, J o. 


(17) 
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When equation (15) is combined with equations (17) and (9a), 


II 

co 2 


~ D °2 C L 

[ L - X .f) 


-Dco 2 c t 

x f 


In 


1 + - 


^°2 C 0 2 

Dco 2 c t 


(18) 


Solving for L/x/- gives 

-l 

(19) 

Next the gas phase diffusivities, A,n 2 > of 0 2 in N 2 and C0 2 in N 2 must be estimated. This is done with the 

Chapman-Enskog correlation [14], which yields diffusivity in cm 2 /s when the variables in the expression below are 
expressed in the appropriate units: 


_L = l+ Dco > C 02 

x f 


Do 2 c t 


In 


1 + 


Dqo 2 c t 


An, =0.0018583 


1 


\V2 f 
7 - 3/2 


2 J 


V 


Pa 2 ., Q 

;,N, ) 


(20) 


Here M-, is the molecular weight (g/mol) of species i, T is the absolute temperature (K), P is the pressure (atm), 
a, jN 2 is the average of the molecular diameter (A) of species i and N 2 , and £\n 2 * s th e collision integral (unitless). 

The later two were determined from the tables of Svehla [15] and Sherwood et al. [16]. Table II lists some 
calculated diffusivities. The value ofi/x/ is also shown in Table II and is constant with temperature. 

For narrow cracks, the contribution of Knudsen diffusion must also be considered. The combined effects of 
ordinary diffusion and Knudsen diffusion have been discussed by other investigators [9-10, 12]. In the Knudsen 
diffusion regime, molecule-wall collisions dominate over molecule-molecule collisions. The effective diffusivity, 
P> eff is a combination of the ordinary diffusivity, D ord , and Knudsen diffusivity, D K . The combined effects of 
ordinary and Knudsen diffusion can be written as [17, 18] 


1 _ 1 1 
/>' rr D old D k 


(21) 


When the mean free path is greater than 10 crack widths, then Knudsen diffusion is the dominating flow mechanism 
[17]. For the 1 atm total pressure used in this study, the mean free path is ~3.3xlCT 7 m. The machined slots are much 
larger than this, and even the craze cracks (-1 ,28x 1 0 5 m) are larger than the mean free path. Thus at the 1 atm (101 
325 Pa) pressure used in this study, the average crack width indicates a very small Knudsen contribution; however, 
at lower pressures which are closer to actual reentry pressure, the Knudsen contribution must be considered. 

Now with values for Xf/L and the gas phase diffusivities carbon consumption can be calculated. As noted, this is 
simply the molar flux of C0 2 (g) emerging from the crack, as given by equations (9b) and (9a): 


J 


1 

co 2 


D C0 2 c co 2 

l - L ~ x f ) 


D 0 2 c 0 2 


1- 


fx f ' 


(22) 


The denominator is written in the above form, since Xf/L is now a known quantity. The goal is to estimate the radius 
of the cavity formed below a crack, as shown in Fig. 6. 

The flux of C0 2 (g) (eq. (11)) can be converted to a mass loss of carbon as 
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M c D o 2 p o 2 wl 


(23) 


dW c 

dt 


M C 

m co 2 



xvl = 




/ M 

M C o 2 RTL 

1- 

V 


Here Wc is the weight loss of carbon, M c is the molecular weight of carbon (12 gm/mol), M^On i s the rnolecular 

weight of carbon dioxide (44 gm/mol), R is the gas constant, w is the crack width, and 1 is the crack length. Note that 
the wl term would be omitted for weight loss per unit area. 

The weight loss with time can be related to the volume change with time of the cavity at the base of the crack. 
Assume a rectangular crack with a half-cylinder cavity of radius r growing at its base. In practice the nonideality of 
the crack is accounted for with a “tortuosity factor” [14]. For this approximation a tortuosity factor of 1 is assumed. 
The weight loss can be related to the radius of the growing half-cylinder cavity as 


d W c dV 
— — = P — 
dt dt 


= P 


nr 2 / 


dt 


(24) 


Here r is the radius of the half cylinder and p is the density of carbon/carbon taken to be 1.362 g/cm 3 [19]. Equating 
equations (24) and (25) and separating variables gives 


) M c Dq 2 Pq 2 wldt 

' M COl PRTL 

1- 

[ Xf \ 
V L )\ 


(25) 


Integrating and solving for r results in the following: 


2 M c Do 2 p 0 2 wt 

M C02 npRTL 

1- 

f Xf ) 


(26) 


Note the crack length cancels out and does not appear in the final expression. Equations (23) and (26) are the key 
equations used to calculate oxidation damage. 

3.2 Diffusion through the Slot and Growing Oxidation Cavity 

Figure 2(a) shows that oxidation through a machined slot has produced a sizable oxidation void. Clearly now 
the model must include transport in the oxidation void as well as the increasing area of reaction. This situation is 
shown schematically in Fig. 6(c). Note also the symmetry of the oxidation void — it is best described with polar 
coordinates. The steady state flux for CCT entering this cylindrical geometry is [20]: 

j;( rJ co 2 h 0 ( 27 ) 

Taking the integration constant to be p, we have: 

r P 

J co, ~~ ( 28 ) 

r 
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We shall simply consider the oxidation cavity as an extension of region II in the slot (Figure 2(a)). Thus we leave 
reaction (4) at the same location in the slot or crack in the SiC coating as in the case for diffusion in the slot or 
crack only. We take the area at the top of cavity which is connected to the crack as A'. It is now more convenient 
to use flux as a weight/(unit time) by multiplying flux in weight/(unit area-unit time) by A'. So in analogy with 

equation (12), we have the following equation for the C0 2 flux entering the oxidation trough, : 


J* A'= -D, 


■CO 2 


co 2 


A' 


r dc co 2 ' 

dr 


A'c„„ J 


CP 2 CO 2 
c T 


(29) 


Equation (28) allows separation of variables: 


A' dr = — A' dr = - D co , A' 


co 2 


( \ 
C T 

C T + c C0 2 


dc 


co 2 


r\ 7"2 

Integrating from r = r, c COl = c QO to r = r 2 c COl = c~ (h gives: 


P^'ln — = -D co A'c T lnl 
r l 1 


C T +c ( 


r 2 

CO 2 


Cf +c 


n 

co 2 


(30) 


(31) 


So the flux at r = r\ is: 



( n ^ 


( r 2 \ 

d C 0 2 A ' c T In 

C T + c C0 2 

D C0 2 nr \ Ic T In 

C T + c co 2 

r\ 

r \ 

P A' 

K c T +c co 2 , 


C T + c co 2 y 



f \ 


f \ 


r 2 


r 2 

rj In 


r\ In 



r l - 

v 1 / 


n 

v 1 y 


(32) 


This is then equated to the expression derived in the previous section for the flux of C0 2 going into region II. We 
are equating an expression in Cartesian coordinates to an expression in polar coordinates and there is a small section 
from x = 0 to r = r\ (see Figure 6(c)), which is ignored in this approximation. The two fluxes, expressions (18) and 
(32), are equated. Note that we have reversed the sign in equation (32) as the polar coordinate system in the 
oxidation void is in the opposite direction from the Cartesian coordinate system in the crack. 


This simplifies to: 


d C 0 2 k1c T In 


D C0 2 c T r l 1 
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n 

C T + C C0 2 

n 


In 
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C T JrC C0 2 
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-nxf 
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' r 2 ^ 

Ct + c co 2 
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c r +c cc>2 J 

r l i 

v 1 y 

n 

, c T* c co 2 , 
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(33) 


In 


(34) 



We need to solve for c^ 0 . In order to do this we need to invert and expand the In function: 


In .y = (x - l)— 2.(x - 1) 2 + ... (2>x>0) 


We use only the first term and obtain: 


f r\ * 

c co 2 ~ c co 2 

In 

' r l' 

_~ nx f 

/'I r 2 

c co 2 - c co 2 

C T + c C0 7 
v y 

, r l. 

n 

, r - 

l CT co 2 J 


Solving for c 


n 

co 2 


gives: 


n 

c co 2 


/ * \ 


f \ 


f ri \ 

c co 2 

In 

r 2 

nx f 

c co 2 
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c T +Cco 2 J 


r \ - 

v 1 / 

n 

n 

{ c t +c co 2 ) 


( \ 


f \ 

1 

In 

r 2 

* 

n 

\ 1 y 

{c T +c C o 2 J 



IX f 


rd cj + c 


n 

co 2 


This expression has the correct limits in that: 


r 2 -» r i 


r i r 2 

c co 2 c co 2 


Xj - — ^ 0 


r\ * 

c co 2 ~^ c co 2 


( 35 ) 


(36) 


(37) 


(38a) 

(38b) 


r\. . ■ r 2 • 

Putting this expression for into the expression for the flux of C0 2 , noting that is close to zero, and 
simplifying gives 


j co 2 a = 


d co 2 kIc co 2 


In 


v 

v 7 " 1 J 


+ (c r + c C q 2 ) 


f n x f ' 

n c r , 


Next we need to determine the change of radius r 2 with time: 


dV __ r 2 d(nr 2 l ) _ M c 1 jtr { _ M c 1 

d co 2 k1c co 2 

dt dt M C q 2 p c °2 M C q 2 p 

In 

f \ 

a 

r \ - 

v 1 y 

+ [ C T +c C0 2 ] 

/ \ 

nxf 

r l C T 
v 1 1 y J 


(39) 


( 40 ) 


Simplifying and separating variables gives: 
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r 2 H>'2 ) - r 2 l n (n ) + r 2 [ C T + c*C0 2 


f nxf^ 
K r l C T JJ 


Me Dcofco^ dt 


M 


co 2 


This can be integrated from t = 0 at r 2 = r\ t = t at r 2 = r 2 , leading the expression: 
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t = - 


co 2 
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c d co 2 c co 2 


2 4 2 2 


X f\c T +cco 2 ) 


r \ c T 


2 2 ( ( * 1^ 

- z ~ z 1 nx f [c 


n__ n_ 

4 2 


*/\ c r +c co 2 
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We can convert this to weight change of the growing half cylinder and obtain the expression: 

w np( r 2~ r lf l 


(41) 


(42) 


(43) 



Equation (42) and (44) replace equations (26) and (24) as expressions for cavity growth when diffusion in both the 
slot or crack and the cavity are considered. 

3.3 Variation of Crack Width due to Thermal Expansion and Contraction 

An important consideration is the variation of the crack width with temperature. This is minor for the machined 
slots, but significant for the craze cracks. The actual coating is not dense SiC, but rather converted carbon-carbon 
fibers and carbon matrix material. For this approximation the accepted coefficient of thermal expansion (CTE) for 
SiC as ~6.1xlCT 6 K -1 is used [21]. The thermal expansion along a carbon fiber is taken as lxlCT 6 K -1 [22]. This puts 
the coating in tension on cool down, which leads to the cracks during processing. Flowever, it leads to compression 
and possible crack closing during heat-up to temperatures above the processing temperature. The extent of crack 
closure can be estimated as 


w RT = Aas(T p-T ) 


(45) 


Flere w RT is the crack width at room temperature, Aa is the difference between the CTE of SiC and that of the 
carbon/carbon, 5 is the crack spacing, T P is the processing temperature (taken as -1650 °C), and T is the operating 
temperature. Since the crack width at room temperature is 1 2.8± 1 .4 pm, the crack spacing is calculated to be 
1.5+0. 2 mm from the equation above, reasonably close to that measured given the approximation of the CTEs. This 
calculated crack spacing is used in the above equation to estimate crack width at the oxidation temperature. At 
1200 °C the crack width was calculated to be 7.1 pm. 

3.4 Effect of Crack Wall Oxidation 

Oxidation along the SiC crack wall may affect the rate of carbon/carbon oxidation. As will be discussed, the 
samples gain considerable weight due to passive oxidation. Microscopy indicates that thin fdms of Si0 2 form on 
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both the internal and external surfaces of the SiC coating. Fig. 7 shows the mouth of a crack and only a thin 
(~0.25-|im) fdm of Si0 2 . 

On the basis of the thin films of oxide observed, the oxidation rate constants for high-purity chemically vapor 
deposited (CVD) SiC [23] are used. The results shown in Table III are somewhat lower than the thicknesses of the 
oxide film on both the slots and cracks. It is probable the lower purity SiC in the conversion coating has an oxidation 
rate constant of an order of magnitude or so greater than that of CVD SiC. Such behavior is well documented [24]. 
Nonetheless images such as those in Fig. 7 indicate that under the conditions of this study, crack closure due to wall 
oxidation can be neglected. 

4. Results and Discussion: Measured Oxidation Kinetics and Comparison with Model 

4.1 Oxidation Below Machined Slots in SiC Coating 

Specimens were oxidized at 1200 °C and examined with both CT and optical microscopy. Cross-sectional 
micrographs for the cavity formed below a 0.37-mm slot is shown in Figs. 2(a). Oxidation damage below all slots 
exhibited this type of half-cylinder morphology. Figs. 8(a) and (b) show the x-ray CT information. 

The data in Figs. 2(a) and 8 make several important points. The similarity between the CT images and the actual 
optical microscopy images give credibility to the use of CT imaging to precisely size and characterize the three- 
dimensional morphology of this type of oxidation damage [25]. Both the optical microscope images and the CT 
images indicated the oxidation damage occurs in the form of a half cylinder. On a macroscopic scale this indicates 
even attack in all directions, supporting diffusion control of the oxidation process and the use of a half cylinder in 
the model. 

Using image analysis software the areas of the oxidation cavities were measured. Approximating these areas as 
semicircles allows calculation of effective radii. The experimental results are compared to the predictions from 
equations (26) and (42) in Figs. 9(a) and (b). Agreement with equation (42) is remarkably good, indicating this is 
the preferred approach to modeling such oxidation behavior. 

A second set of experiments was done to follow the kinetics of oxidation below machined slots with weight loss 
as the indicator of oxidation damage. Kinetics were followed using interrupted weight loss measurements. A typical 
kinetic curve is shown in Fig. 10. Note there is an initial weight gain, which is likely passive oxidation of the SiC, 
followed by a weight loss. The passive oxidation of the SiC will be discussed in the following section. The rates of 
weight loss are the carbon oxidation rates per the exposed area of carbon, which was determined from the slot 
dimensions. Note that the experimental weight loss rates appear to be roughly linear. It was found that the linear 
terms in equation (44) tended to dominate and thus the results from equation (44) were plotted and fitted to a linear 
rate for comparison to experiment. 

Table IV compares the measured rates with those predicted from equation (44). Rates are normalized per unit 
area of exposed carbon, and hence the predictions are all the same for the given temperature of 1200 °C. The error in 
the predicted oxidation rate values is due to uncertainty in the measurement of slot depth (coating thickness). The 
calculated rates are within a factor of 1 to 2.6x the measured, which is reasonable. 

4.2 Oxidation Below Craze Cracks in the SiC Coating 

Having established the validity of the model with the machined slots, it is now appropriate to proceed to the 
more realistic situation of craze cracks in the SiC coating. Now there is a less well-defined geometry, as illustrated 
with the crack pattern in Fig. 4 and crack width measurement in Fig. 5. Nonetheless these measurements can be used 
to estimate approximate values for the diffusion model. 

Oxidation kinetics for these specimens were followed with a recording thermobalance that enabled continuous 
weight change measurements. Results for 1000, 1 100, 1200, and 1300 °C are given in Figs. 1 1(a) to (d). These data 
are reported in weight change per geometrical surface area. As will be discussed, the geometric surface area must be 
modified for actual reactive area of oxidation, which is different for each of the two processes occurring. 

The oxidation kinetics suggest a combined parabolic/linear rate law. The initial weight gain observed at the 
three higher temperatures is likely due to the passive oxidation of SiC, which is parabolic: 
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is the specific weight change (weight per unit area), k p is the parabolic rate constant, and t is the time. 

The specific weight change over the surface area A of the SiC is taken as the geometric area of the sample disk. The 
weight loss after ~0.5 hr is linear with time and most likely due to carbon oxidation through the craze cracks. It is 
described by 

AM'] , 

— J = ¥ ( 4? ) 

Here k;is the linear rate constant. Note that A' is the exposed area of carbon/carbon. This is the crack width at 
temperature multiplied by the total crack length for a particular sample. 

The changes in the kinetic curves with temperature support the interpretation of the process as passive SiC 
oxidation coupled with C/C oxidation. At 1000 °C the initial parabolic regime was not observed; however, at 1200 
and 1300 °C it was quite evident. This is due to the fact that passive oxidation of SiC shows a strong dependence on 
temperature [23], Note, however, that the linear rate constant is relatively temperature independent. Diffusion- 
controlled oxidation of carbon shows only a weak temperature dependence, a combination of the T 1 ' 2 dependence of 
the diffusion coefficient and the T 1 dependence of the crack closure due to thermal expansion. This dependence may 
not be evident in the limited data set in this study. 

The combined passive oxidation of SiC/linear oxidation of carbon rate law is written as 

[^f\ =Vv+^ ( 48 ) 

v A ) total 



Here 



is the measured weight change per unit of geometrical area, k p is the parabolic rate constant, t is 

total 


the time, and kj is the linear rate constant. This expression is similar to that used to describe the combined 
oxidation/vaporization of chromium alloys [26] and SiC [27], but this is somewhat simpler in that there are two 
independent processes occurring. 

The dashed lines in Figs. 1 1(a) to (d) are fits of equation (48) to each kinetic curve. From these fits the parabolic 
and linear rate constants can be derived. Fits to the linear region are good; fits to the parabolic region are only 
approximate because of such factors as changing surface areas of passive oxidation and passive oxidation internal to 
the coating, which are not included in equation (48). The parabolic and linear rate constants extracted from the fit to 
the kinetic data are given in Table V. The parabolic rate constant was taken directly from the fit. Note that the listed 
parabolic rate constants are many orders of magnitude greater than those measured for high-purity SiC [23]. This is 
due to extensive internal oxidation of the porous SiC conversion coating. Measurements of the SiC conversion 
coating density have been reported to be 2.1 14 g/cm 3 [19]. When compared with the theoretical density of 3.217 
g/cm 3 for SiC, this indicates a high amount of porosity. This porosity leads to the large extent of internal oxidation. 

The linear rate constant shown in Table V was taken from the fit and multiplied by the geometrical surface area. 
This gives a rate of carbon oxidation per area of exposed carbon/carbon and allows direct comparison to the 
calculated rate in Table V. The predicted rates were obtained from equation (24). 

The cracks in the SiC form a complex geometry. The oxidation cavities are an irregular shape as seen in the 
micrograph in Fig. 2(b). Fig. 12(a) shows the location of two CT scans, and Fig. 12(b) shows the results of the CT 
scans illustrating the varying size and shape of the oxidation cavity below a craze crack. The actual oxidation 
cavities are darker and are outlined in white in Fig. 12(b). 

Despite this complex geometry of both the craze crack and the resultant oxidation damage, average values of 
crack width and length can be put into the oxidation model to estimate the rates. These are shown in Table V. The 
predicted and measured rates show very good agreement at 1000 and 1100 °C, which is remarkable given the 
irregular geometry of the cracks and corresponding oxidation damage as well as the approximations used in this 
study. We believe the larger variation at 1200 and 1300 °C may be due to volatile species in the conversion coating 
which the net weight loss measurement would necessarily include [4]. Observation of the both the 1000 and 1300 
°C specimens after 2.5 hr indicate that the oxidation voids below the cracks are no longer distinct but start to merge, 
indicating that the approximation of the oxidation process as a single trough is no longer valid. This may be another 
reason for the discrepancy between the calculated and measured rates. 
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5. Conclusions 


The kinetics of subsurface oxidation of SiC-protected reinforced carbon/carbon has been studied in the diffusion 
control regime (1000 to 1300 °C) in air. Specimens with machined slots in the SiC coating as well as specimens 
with natural occurring craze cracks in the SiC coating were tested. The geometry of these oxidation pathways was 
characterized with the aim of obtaining the exposed area of the carbon/carbon. For the craze cracks, crack closure 
due to shrinkage was considered. Oxidation of the SiC walls was shown to be negligible under the conditions of this 
study. Oxidation damage of the carbon/carbon was assessed with weight changes, CT scans, and optical microscopy 
of sections. The CT scans confirm the idealized half-cylinder shape of the oxidation cavity in the case of a slot. A 
two-step diffusion control model for carbon oxidation is adapted to describe the growth kinetics of this half cylinder. 
Two models were used — a simplified model which considered diffusion only in the slot or craze crack and a more 
detailed model which considered diffusion in both the slot or craze crack and the growing oxidation cavity. 

Oxidation rates were first measured from optical microscopy of the growing oxidation cavity for the well- 
defined geometry of the slotted specimens. For these optical measurements, the data showed remarkably good 
agreement with the more detailed model which considered diffusion in the slot and the oxidation cavity. A series of 
experiments were also done where oxidation rates were followed with overall weight loss for both the specimens 
with slots and natural crazed cracks. For the weight loss measurements agreement to the model which considered 
diffusion in both the slot or craze crack and the growing oxidation cavity was reasonable. 
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Appendix — Symbols 


A 

A’ 

c 

D 

J 

k 

L 

l 

M 

P 

q 

R 

r 

s 

T 

t 

v 

W 

w 

X 

x f 

CT 

Q 

a 

P 

Subscripts: 

/ 

1 
1 

P 

P 

T 

1 

2 

Superscripts: 

ave 

eff 

K 

L 

ord 

0 

I 

II 
* 


area of SiC for oxidation 

area of exposed carbon-carbon for oxidation 

molar concentration 

diffusivity 

molar flux 

rate constant (weight 2 /(unit area) 2 -time) 
depth of crack or coating thickness 
crack length 

molecular weight (g/mol) 
total pressure in atmospheres 
pore radius 
gas constant 

radius of half-cylinder cavity 
crack spacing 
temperature 
time 

molar velocity 
weight loss 
crack width 

depth of crack coordinate with origin at SiC/C-C interface 
depth to C0/C0 2 (region II) and C0 2 /0 2 (region I) boundary 
average collision diameter between species i and N 2 (angstroms) 
collision integral 
coefficient of thermal expansion 

integration constant for steady flux in cylindrical symmetry 

indicates boundary of regions I and II 

species 

linear 

processing 

parabolic 

total 

designates radius of half cylinder cavity before oxidation (taken as 0.5*crack width) 
designates radius of half cylinder cavity after oxidation 

average 

effective 

Knudsen 

indicates position L 
ordinary 

indicates position 0 
region 1 
region II 

indicates position Xf 
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Figure Captions 

Fig. 1. — Schematic of reinforced carbon/carbon (RCC). 

Fig. 2. — Optical micrographs showing oxidation void in RCC formed below crack in SiC coating for (a) 0.37-mm 
slot in air at 1200 °C for 2.5 hr (b) 1 100 °C and 0.0066 atm (667 Pa) for 1 hr. 

Fig. 3. — Machined slots in RCC buttons. 

Fig. 4. — Crack pattern on RCC disk, (a) Optical micrograph, (b) “Skeleton” trace of cracks. 

Fig. 5. — Crack in SiC conversion coating. 

Fig. 6. — Two-step diffusion control model, (a) Schematic of fluxes and boundaries for a slot or crack, (b) Relative 
concentrations of gases, (c) Schematic of idealized void. 

Fig. 7. — Oxide fdm on surface of SiC coating (0.5 hr in air at 1200 °C) at crack opening, (a) Electron micrograph, 
(b) Oxygen map. 

Fig. 8. — Computed tomography (CT) of machine-slotted SiC-coated RCC oxidized for 0.5 hr in air at 1200 °C. (a) 
CT sample locations, (b) Three CT slices. Slot near notch on left is 0.76 mm, and slot on right is 1.03 mm. (c) 
Three-dimensional reconstruction of disk. 

Fig. 9. — Oxide cavity growth, (a) 0.53 mm slot, (b) 1.05 mm slot. Both plots show the calculated results from 
considering diffusion only in the slot and diffusion on the slot and growing cavity. 

Fig. 10. — Weight loss oxidation kinetics determined from interrupted measurements for RCC specimen with 0.484- 
mm slot width at 1200 °C in air. 

Fig. 11.— Weight loss oxidation kinetics for RCC in air. (a) 1000 °C. (b) 1100 °C. (c) 1200 °C. (d) 1300 °C. 

Fig. 12. — Computed tomography (CT) of craze crack of SiC-coated RCC oxidized for 2.5 hr in air at 1143 °C. (a) 
Location of CT slices, (b) Two CT slices. 
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TABLE I. Measured crack parameters for SiC-coated carbon/carbon materials 


SiC-coated carbon/carbon 
material 

Coating 

thickness, 

mm 

Crack length per 
unit area, 
mm 1 

Crack spacing, 
mm 

Crack width, 
pm 

Crack area per 
unit area 

Present study, RCC 
Yurgartis et al. [2], Specimen A 
Yurgartis et al. [2], Specimen B 
Buchanan and Little [3] 

0.78+0.14 

0.3024+0.03866 

0.2788+0.02853 

0.33+0.04 

2.985 

3.555 

3. 4+0.9 

0.6247+0.2233“ 

0.6517+0.3356“ 

12.8+1.41 

1.72+1.41 

2.02+1.47 

4.2x1 0‘ 3 
5.1xl0 -3 
7.3xl0“ 3 
(3-8)xl0^ 


a Through coating thickness. 
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TABLE II. Diffusivities of 0 2 (g) and C0 2 (g) in N 2 (g) and 
position of boundary between regions I and II (Fig. 6(a)). 
[Total pressure is 1 atm (101 325 Pa).] 


Temperature, 

°C 

D 0„N 2 > 

cm 2 /s 

A:o 2 ,n 2 , 
cm 2 /s 

L/jCf 

1000 

2.38 

1.86 

2.13 

1100 

2.70 

2.11 

2.13 

1200 

3.03 

2.38 

2.13 

1300 

3.38 

2.65 

2.13 
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TABLE III. Estimated oxide film thicknesses on slot or crack walls of SiC coating 


Temperature, 

°C 

Rate constant 

Oxide thickness | 

1 atm 0 2 , a 
pm 2 /hr 

0.21 atm 0 2 , 
pm 2 /hr 

0.5 hr in 0.21 atm 0 2 , 
pm 

2.5 hr in 0.21 atm 0 2 , 
pm 

1200 

1.82xl0“ 2 

3.78xl0' 3 

4.35xl0~ 2 

9.72xl0“ 2 

1300 

3.76xl0“ 2 

7.90x1 O' 3 

6.28xl0 -2 

0.14 


a From reference 23. 
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TABLE IV. Oxidation rates for SiC-coated 
RCC specimens with machined slots 

[At T= 1200 °C in air.] 


Slot width, 
w, 
mm 

Slot length, 
/, 

mm 

Weight loss rate, 
mg/mm 2 -hr 

Measured 

Calculated 3 

0.484+0.062 

8.177+0.923 

34+6 

18.3 

0.312+0.033 

8.284+0.63 

48+6 

20.4 

0.560+0.013 

7.094+0.132 

48+1 

17.6 

0.466+0.04 

7.114+0.019 

57+5 

18.6 

1.116+0.024 

7.602+0.101 

18+1 

13.9 

0.963+0.0041 

7.624+0.039 

21+1 

14.7 


“From equation (44). 
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TABLE V. Measured and calculated rates of RCC oxidation through craze cracks 


Sample 

temperature, 

°C 

Geometrical 
surface area, 
mm 2 

Total crack 
length, 3 
mm 

Crack 

width, 

pm 

Area of carbon 
exposed by 
craze cracks, b 
mm 2 

Calculated 
linear rate, c 
mg/mm 2 -hr 

Measured 
linear rate, d 
mg/mm 2 -hr 

Measured 
parabolic rate, 6 
mg 2 /mm 4 -hr 

Room temperature 



12.8+1.4 





1000 

850.2+10 

281+34 

10.2+2.7 

2.87+0.8 

25 

31+9 


1100 

935.5±10 

309+37 

8. 7+2. 3 

2.67+0.8 

23.3 

31+9 

3.94xl0“ 3 

1200 

946.8±10 

312+37 

7.1+1. 8 

2.21+0.6 

21.8 

53+15 

2.79xl0“ 2 

1300 

940.2+10 

310+37 

5. 5+1. 4 

1.71+0.5 

20.5 

67+20 

3.25xl0“ 2 


a Equal to (geometrical area x crack length)/unit area. 

b Equal to (total crack length) x (crack width). 

c From equation (44). 

d Area is exposed area of carbon/carbon. 

c Area is geometrical surface area. 
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